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Abstract 

We investigate the decay of a scalar field outside a Schwarzschild anti de Sitter black 
hole. This is determined by computing the complex frequencies associated with quasinor- 
mal modes. There are qualitative differences from the asymptotically flat case, even in 
the limit of small black holes. In particular, for a given angular dependence, the decay is 
always exponential - there are no power law tails at late times. In terms of the AdS/CFT 
correspondence, a large black hole corresponds to an approximately thermal state in the 
field theory, and the decay of the scalar field corresponds to the decay of a perturbation 
of this state. Thus one obtains the timescale for the approach to thermal equilibrium. 
We compute these timescales for the strongly coupled field theories in three, four, and six 
dimensions which are dual to string theory in asymptotically AdS spacetimes. 
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1. Introduction 



It is well known that if you perturb a black hole, the surrounding geometry will "ring" , 
i.e., undergo damped oscillations. The frequencies and damping times of these oscillations 
are entirely fixed by the black hole, and are independent of the initial perturbation. These 
oscillations are similar to normal modes of a closed system. However, since the field can 
fall into the black hole or radiate to infinity, the modes decay and the corresponding 
frequencies are complex. These oscillations are known as "quasinormal modes" . For black 
holes in asymptotically flat spacetimes, they have been studied for almost thirty years 
0,3- The radiation associated with these modes is expected to be seen with gravitational 
wave detectors in the coming decade. Motivated by inflation, the quasinormal modes of 
black holes in de Sitter space have recently been studied [[|,f|] . 

For spacetimes which asymptotically approach anti de Sitter (AdS) spacetime the 
situation is slightly different. In the absence of a black hole, most fields propagating in 
AdS can be expanded in ordinary normal modes. The cosmological constant provides an 
effective confining box, and solutions only exist with discrete (real) frequencies. However 
once a black hole is present, this is no longer the case. The fields can now fall into the 
black hole and decay. There should exist complex frequencies, characteristic of the black 
hole, which describe the decay of perturbations outside the horizon. We will compute these 
quasinormal frequencies below, for spacetimes of various dimensions. 

The quasinormal frequencies of AdS black holes have a direct interpretation in terms 
of the dual conformal field theory (CFT) According to the AdS/CFT corre- 

spondence, a large static black hole in AdS corresponds to an (approximately) thermal 
state in the CFT. Perturbing the black hole corresponds to perturbing this thermal state, 
and the decay of the perturbation describes the return to thermal equilibrium. So we 
obtain a prediction for the thermalization timescale in the strongly coupled CFT. It seems 
difficult to compute this timescale directly in the CFT. Since the system will clearly not 
thermalize in the free field limit, at weak coupling, this timescale will be very long and 
depend on the coupling constant. In the limit of strong coupling, it seems plausible that 
the timescale will remain nonzero, and be independent of the coupling. This is because 



1 The importance of these modes in AdS was independently recognized in B, but they were 



not computed. They were computed in [1C], but only for a conformally invariant scalar field whose 
asymptotic behavior is similar to flat spacetime. The confining behavior of AdS is crucial for the 
AdS/CFT correspondence. 
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the initial state is characterized by excitations with size of order the thermal wavelength, 
so causality suggests that the relaxation timescale should also be of order the thermal 
wavelength. 

The results we obtain are consistent with this expectation. A black hole in AdS is 
determined by two dimensionful parameters, the AdS radius R and the black hole radius 
r + . The quasinormal frequencies must be functions of these parameters. For large black 
holes, r + ^> R, we will show that there is an additional symmetry which insures that 
the frequencies can depend only on the black hole temperature T ~ r + /R 2 . However, 
for smaller black holes, this is no longer the case. Whereas the temperature begins to 
increase as one decreases r_|_ below R, we find that the (imaginary part of the) frequency 
continues to decrease with r+. This is different from what happens for asymptotically flat 
black holes. An ordinary Schwarzschild black hole has only one dimensionful parameter 
which can be taken to be the temperature. Its quasinormal frequencies must therefore be 
multiples of this temperature. Thus small black holes in AdS do NOT behave like black 
holes in asymptotically flat spacetime. The reason is simply that the boundary conditions 
at infinity are changed. More physically, the late time behavior of the field is affected by 
waves bouncing off the potential at large r. 

Another difference from the asymptotically flat case concerns the decay at very late 
times. For a Schwarzschild black hole, it is known that the exponential decay associated 
with the quasinormal modes eventually gives way to a power law tail This has been 
shown to be associated with the scattering of the field off the Coulomb potential at large 
r. As we will discuss later, for asymptotically AdS black holes, this does not occur. 

We will compute the quasinormal frequencies for Schwarzschild-AdS black holes in 
the dimensions of interest for the AdS/CFT correspondence: four, five, and seven. We will 
consider minimally coupled scalar perturbations representing, e.g., the dilaton. This cor- 
responds to a particular perturbation of the CFT. For example, for AdS$, it corresponds 
to a perturbation of an (approximately) thermal state in super Yang-Mills on S 3 x R with 
< F 2 > nonzero. In the linearized approximation we are using, the spacetime metric is 
not affected by the scalar field. So the perturbation of the thermal state does not change 
the energy density, which remains uniform over the sphere. The late time decay of this 
perturbation is universal in the sense that all solutions for the dilaton with the same angu- 
lar dependence will decay at the same rate, which is determined by the imaginary part of 
the lowest quasinormal frequency. Different perturbations, corresponding to different lin- 
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earized supergravity fields, will have different quasinormal frequencies and hence decay at 
different rates. Although we work in the classical supergravity limit, our results would not 
be affected if one includes small semiclassical corrections such as black holes in equilibrium 
with their Hawking radiation. 

A brief outline of this paper is the following. In the next section we review the 
definition of quasinormal modes, their relation to the late time behavior of the field, and 
derive some of their properties using analytic arguments. The numerical approach we use 
to compute the complex frequencies is described in section 3. In the following section we 
discuss the results for both large black holes r + ^> R and intermediate size black holes 
r + ~ R. In section 5 we consider the limit of small black holes r + <C R. Although there 
is a striking similarity between some of our results and some results obtained in the study 
of black hole critical phenomena [[12], we will argue that this is probably just a numerical 
coincidence. The conclusion contains some speculations about the CFT interpretation of 
the quasinormal frequencies in the regime where they do not scale with the temperature. 
In the appendix, we give some more details on our numerical calculations. 

2. Definition of quasinormal modes and analytic arguments 

Since we are interested in studying AdS black holes in various dimensions, we begin 
with the d dimensional Schwarzschild-AdS metric: 





where 




(2.2) 



R is the AdS radius and ro is related to the black hole mass via 



(d - 2)A d - 2 r, 
167rG d 



d-3 



M = 



(2.3) 



where Aa-2 = 2% 2 /r(i_i) j s the area of a unit (d — 2)-sphere. The black hole horizon 
is at r = r+, the largest zero of /, and its Hawking temperature is 



T 



f(r+) _ (d-l)rl + (d-3)R 2 
47r 4irr + R 2 



(2.4) 
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We are interested in solutions to the minimally coupled scalar wave equation 



V 2 $ = (2.5) 

If we consider modes 

<&(£, r, angles) = r 1 ^ tp(r) Y(angles) e~ iujt (2.6) 

where Y denotes the spherical harmonics on S d ~ 2 , and introduce a new radial coordinate 
dr* = dr I f(r), the wave equation reduces to the standard form 

{dl+u 2 -V(r*)]i; = 0. (2.7) 

The potential V is positive and vanishes at the horizon, which corresponds to r* = — oo. 
It diverges at r = oo, which corresponds to a finite value of r*. 

To define quasinormal modes, let us first consider the case of a simple Schwarzschild 
black hole. Since the spacetime is asymptotically flat, the potential now vanishes near 
infinity. Clearly, a solution exists for each uj corresponding to a wave coming in from 
infinity, scattering off the potential and being partly reflected and partly absorbed by the 
black hole. Quasinormal modes are defined as solutions which are purely outgoing near 
infinity $ ~ e — iw(t— r«) anc j p Ure iy ingoing near the horizon $ ~ e -ioj(t+r*) j\j j n jtj a i 

incoming wave from infinity is allowed. This will only be possible for a discrete set of 
complex uj called the quasinormal frequencies. 

For the asymptotically AdS case, the potential diverges at infinity, so we must require 
that $ vanish there. In the absence of a black hole, r* has only a finite range and solutions 
exist for only a discrete set of real uj. However once the black hole is added, there are 
again solutions with any value of uj. These correspond to an outgoing wave coming from 
the (past) horizon, scattering off the potential and becoming an ingoing wave entering the 
(future) horizon. Quasinormal modes are defined to be modes with only ingoing waves 
near the horizon. These again exist for only a discrete set of complex uj. 

It should perhaps be emphasized that these modes are not the same as the ones that 



have recently been computed in connection with the glueball masses |T3[ . There are several 
differences: First, the background for the glueball mass calculation is not the spherically 
symmetric AdS black hole, but an analytic continuation of the plane symmetric AdS black 
hole. Second, because of the analytic continuation, the horizon becomes a regular origin, 
and the boundary conditions there are not the analytic continuation of the ingoing wave 
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boundary condition imposed for quasinormal modes. Finally, the glueball masses are real 
quantities, while as we have said, the quasinormal frequencies will be complex. This makes 
them more difficult to compute numerically. 

One can show that the complex quasinormal frequencies determine the fall off 

of the field at late times. The basic idea is to start by writing the solution to the wave 
equation in terms of the retarded Green's function and initial data on a constant t surface. 
One then rewrites the Green's function in terms of its Fourier transform with respect to 
t. The quasinormal modes arise as poles of the Green's function in the complex frequency 
plane, and their contributions to the solution can be extracted by closing the contour with 
a large semicircle near infinity. 

For a black hole in asymptotically flat spacetimes, Price []T]J showed that after the 
exponential decay due to the quasinormal ringing, the field will decay as a power law 
$ ~ f-(2i+3) w h ere i i s th e angular quantum number. This has been seen explicitly in 



numerical simulations [[15). Mathematically, this is due to a cut in the Green's function 
along the negative imaginary frequency axis. More physically, this behavior is due to 
scattering off the weak Coulomb potential near infinity. For the case of a black hole in 
AdS, the potential diverges at infinity and vanishes exponentially near the horizon. Ching 
et. al. |14| have analyzed the late time behavior of a broad class of wave equations with 
potentials. They show that there are no power law tails for a potential which vanishes 
exponentially. So there will be no power law tails for black holes in AdS. 

For a black hole with radius much smaller than the AdS radius, one might expect 
an intermediate time regime where one sees power law behavior before the new boundary 
conditions at infinity become important. However, this would occur only if one starts 
with large quasinormal modes with 10 ~ 1/V+ associated with a Schwarzschild black hole. 
We will see that the lowest modes of a Schwarzschild-AdS black hole are much smaller 
and their exponential decay is so slow that it eliminates the intermediate time power law 
behavior. 

The quasinormal frequencies will in general depend on the two parameters in the 
problem R, tq. By rescaling the metric, ds = A 2 (is 2 , and rescaling the coordinates i = Xt 
and f = Ar, the new metric again takes the form ( |2.1|) with rescaled constants R and 
tq. Since the wave equation ( |2.5| ) is clearly invariant under this constant rescaling of the 
metric, we can use it to set e.g. R = 1. This rescaling is possible for any metric and 
physically just corresponds to a choice of units. In our case, we measure all quantities in 
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units of the AdS radius. The quasinormal frequencies can still be arbitrary functions of 

ro- 

We now show that for large black holes, ro ^> R, the frequencies must be proportional 
to the black hole temperature. This is a result of an independent scaling one can do in 
this limit. For large black holes, the region outside the horizon of the Schwarzschild-AdS 
metric ( |2.1| ) becomes approximately plane symmetric: 

ds 2 = -h(r) dt 2 + fc(r) _1 dr 2 + r 2 dx t dx l (2.8) 

where 

%) = 4-(-V" 3 - (2-9) 



R? \ r 

For this metric one can rescale ro by a pure coordinate transformation: t = at, x% = 
a&i, r = f/a for constant a. This does not rescale the overall metric, or the AdS radius 
R. The horizon radius r^T 1 = R 2 r^~ 3 gets rescaled by r + = r+/a. Of course, under this 
coordinate transformation of the metric, solutions of the wave equation are related by the 
same coordinate transformation. For solutions which are independent of x l (the analog 
of the I = modes) we have e _lw ( r+ ' >i = e _ltt;( - r+ - )t , which implies u(r+) oc r+. Since the 
Hawking temperature of the metric ( |2.8|) is also proportional to the horizon radius, 

m d — 1 r_i_ . 

< 2 - 10) 

we see that the frequencies must scale with the temperature for large black holes. For 
solutions proportional to e lkiX , this scaling argument implies u(ar+,aki) = au)(r+,ki). 
So if ^> kik 1 , one can rescale so that k 2 is negligibly small. The above argument then 
shows that ui still scales with the temperature. One can then rescale back to r+ > J? 
to apply to large black holes. In other words, for any fcj, the quasinormal frequencies 
scale with the temperature in the limit of large temperatures T 2 ^> k 2 . This argument 
does not apply to black holes of order the AdS radius, and indeed we will find that the 
quasinormal frequencies do not scale with the temperature in this regime. But it does 
confirm the expectation that the approach to thermal equilibrium in the dual field theory 
should depend only on the temperature (at least for large temperature). 

Since we want modes which behave like e - tU} ( t + r *) near the horizon, it is convenient to 
set v = t+r*, and work with ingoing Eddington coordinates. The metric for Schwarzschild- 
AdS in d dimensions in ingoing Eddington coordinates is 

ds 2 = -f(r)dv 2 + 2dvdr + r 2 dfl 2 d _ 2 (2.11) 
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where / is again given by (|2.2| ) . The minimally-coupled scalar wave equation ( |2.5| ) may be 
reduced to an ordinary, second order, linear differential equation in r by the separation of 
variables, 

r, angles) = r^ ip(r) Y(angles) e~ iujv (2.12) 
This yields the following radial equation for ip(r): 

/(r ) ±_^ r ) + [f (r) - 2iu] ^-VW - F(r) ^(r) = 0, (2.13) 

with the effective potential V(r) given by (R = 1) 

(d-2)(d-4) „ ,,d-2 c 
V ( r ) = + ~^rf ( r ) + ^ 

_ d(d-2) (d-2)(d-4)+4c (ci-2) 2 r rf - 3 

where 

c = Z(Z + - 3) (2.15) 

is the eigenvalue of the Laplacian on S d ~ 2 . Note that V(r) is manifestly positive for d > 4. 

Ingoing modes near the (future) horizon are described, of course, by a nonzero multiple 
of e~ lLJV . Outgoing modes near the horizon can also be expressed in terms of ingoing 
Eddington coordinates via e~ lU}<Kt ~ r ^ = e - tuJV e 2lU}r * . Since 

dv 1 

ln(r-r+) (2.16) 



f(r) f'(r+) 

near the horizon r = r + , the outgoing modes behave like 

;(t-r„) _ e —iuv e 2iiur* ^ g-iwv^. _ \2iu/ f'(r+) (2.17) 



Since v, r are good coordinates near the horizon, the outgoing modes are not smooth (C°°) 
at r = r_|_ unless 2iu)/ f (r+) is a positive integer. We show below that the imaginary part 
of u must be negative, so the exponent in ( |2.17|) always has a positive real part. Thus 
the outgoing modes vanish near the future horizon, while the ingoing modes are nonzero 
there. However we also show (in the next section) that 2iuj / f'(r + ) cannot be a positive 
integer, so the outgoing modes are not smooth at r = r + . 

We wish to find the complex values of uj such that ( |2.13| ) has a solution with only 
ingoing modes near the horizon, and vanishing at infinity. We will eliminate the outgoing 
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modes by first assuming the solution is smooth at r = r + , and then showing that the 
allowed descrete values of u are such that 2iu/ f'(r + ) is not an integer. The actual values 
of 10 must be computed numerically, but some general properties can be seen analytically. 
For example, we now show that there are no solutions with iui pure real, and 2iu < f'(r+). 
If iu were real, then the equation would be real and the solutions ip would be real. If there 
were a local extremum at some point f, then ip'(r) = and ip"(f) would have the same 
sign as i/j(r). So if ip were positive at r, it would have to increase as r increased. Similarly, 
if it were negative, it would have to decrease. In neither case, could it approach zero 
asymptotically. We conclude that the solutions must monotonically approach zero. Now if 
2iu> < /'(r+), i/j'(r+) has the same sign as ?/>(r_|_).ll So as one moves away from the horizon, 
the solutions move farther away from zero and hence can never reach zero asymptotically. 
This analytic argument only applies if 2iuj < f'(r + ). But we will see numerically that even 
without this restriction, there are no solutions with iu pure real. 

A more powerful result can be obtained as follows. Multiplying ( |2.13| ) by ip and 
integrating from r_|_ to oo yields 



dr 

r+ 



dr \ dr I dr 



(2.18) 



The first term can be integrated by parts without picking up a surface term since f(r+) = 
and ^(oo) = 0. This yields 



/•OO 

/ dr[f\ij'\ 2 + 2iu^' + V\ip\ 2 } = (2.19) 

Jr + 

Taking the imaginary part of ( |2.19| ) yields 

dr[u$ip' +unp$'\ = (2.20) 



Integrating the second term by parts yields 



{u-Q) I drW = w|^(r + )| 2 (2.21) 

Jr+ 



2 This is where the condition of no outgoing modes near the horizon is used. If outgoing waves 
were present, f(r)ip"(r) would no longer vanish at r = r+, and V/( r +) need not have the same 
sign as tp(r+). 
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Substituting this back into ( |2. 19|) we obtain the final result 

jT^' + ^=- M (2.22) 

Since / and V are both positive definite outside the horizon, this equation clearly shows 
that there are no solutions with Im 10 > 0. These would correspond to unstable modes 
which grow exponentially in time. There are also no solutions with Im uj = 0: All solutions 
must decay in time. In addition, eq. ( |2.22| ) shows that the only solution which vanishes at 
the horizon (and infinity) is zero everywhere. Since the equation is linear, we can always 
rescale ijj so that i/i(r+) = 1. 



3. Numerical approach to computing quasinormal modes 

To compute the quasinormal modes, we will expand the solution in a power series 
about the horizon and impose the boundary condition that the solution vanish at infinity. 
In order to map the entire region of interest, r+ < r < oo, into a finite parameter range, 
we change variables to x = 1/r. In general, a power series expansion will have a radius 
of convergence at least as large as the distance to the nearest pole. Examining the pole 
structure of ( [2.13| ) in the whole complex r plane, we find d + 1 regular singular points, 
at r = 0, r = oo, and at the d — 1 zeros of /, one of which, r = r + corresponds to the 
horizon. At least for d = 4, 5 or 7, if we use the variable x = 1/r and expand about the 
horizon, x+ = l/r+, the radius of convergence will! reach to x = 0, so that we can use 
this expansion to consider the behavior of the solution as r — > oo. 

In terms of our new variable x = 1/r, ( |2.13| ) becomes 

s(x) ^(x) + ^(x) + U(X) iP(x) = (3.1) 

dx z x — X+ dx [x — x+y 

where the coefficient functions are given by 

d-3 d+l 4 _ 2 r 2 , -1 2 _i_ i i i 

s(x) = — = + . = x d + . . . + — S — x 4 + -^-x 6 H x 2 (3.2) 

x — x + x + x% x+ x+ 

t(x) = (d - 1) r^" 3 x d - 2x 3 - 2x 2 iu (3.3) 



3 For d = 4 and d = 5, one can show analytically that starting at the horizon, x = x+, the 
nearest pole is indeed x = 0. For d = 7 we have checked numerically that this is again the case. 
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u(x) = (x-x+)V(x) (3.4) 
The parameter r ~ should be viewed as a function of the horizon radius: r ~ — — 



d-l 



Since s, t, and u are all polynomials of degree d we may expand them about the horizon 
x = s(x) = Ylt,=o Sn ( x ~ x +) n i an d similarly for t(x) and u(x). It will be useful to 
note that sq = 2x^_ k, to = 2x^_ (k — iu), and uq = 0, where k is the surface gravity, which 
is related to the black hole temperature (|2.4|) by 



K = = 2nT. (3.5) 



Also, since so ^ 0, x = x + is a regular singular point of (|3TT|). 

To determine the behavior of the solutions near the horizon, we first set ip{x) = 



(x — £-|_) Q and substitute into (|3.1|) . Then to leading order we get 



a(a — 1) Sq + a to = 2x^_a («k- iu) = (3-6) 



which has two solutions a = and a = iu/n. We see from ( 2.17|) that these correspond 



precisely to the ingoing and outgoing modes near the horizon respectively. Since we want 
to include only the ingoing modes, we take a = 0. This corresponds to looking for a 
solution of the form 

oo 

^(x) = ^a n (x-x+) n (3.7) 



n=0 



Substituting ( |3.7| ) into ( |3.1| ) and equating coefficients of (x — x + ) n for each n, we obtain 
the following recursion relationsS for the a n : 

^ n — 1 

a n = -— 2J [^(fe - 1) S n-fc + ^n-fc + Un-k] Clk (3.8) 
^ n fc=0 

where 

P n = n(n — 1) so + n t = 2x+n (n k — iuj) (3-9) 

Since the leading coefficient ao is undetermined, this yields a one parameter family of 
solutions, as expected for a linear equation. 



4 Although the standard way of writing ( |3.1| ) is to set the coefficient of ip" to 1 which yields 
simpler-looking recursion relations, the advantage of the present formulation is that since s(x), 
t(x), and u(x) are polynomials, their analytic expansions will terminate after a finite number of 
terms, so that each a n will be given in terms of a relatively small number of terms. 
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The solutions to (|2.5| ) in asymptotically AdS spacetime are $ ~ constant and <E> ~ 
\jr d ~ x as r — > oo, which translates into ^ ~ and ?/> ~ r _d / 2 , respectively. We are 
interested in normalizable modes, so we must select only solutions which satisfy ijj — > as 
r — > oo (or x — > 0). This means that we require (|3.7| ) to vanish at x = 0, which is satisfied 
only for special (discrete) values of u. (For all other values of u, the solution will blow 
up, if)(0) = oo.) Thus in order to find the quasinormal modes, we need to find the zeros of 
^2^=0 a n(^) (— x+) n in the complex u plane. This is done by truncating the series after a 
large number of terms and computing the partial sum as a function of uj. One can then find 
zeros of this partial sum, and check the accuracy by seeing how much the location of the 
zero changes as one goes to higher partial sums. Some details are given in the Appendix. 

One can now easily show that 2iw// / (r_|_) = iu/n cannot be an integer. If ui is 
pure imaginary and iui = fin for some integer n, then = 0. This implies an additional 
constraint on the coefficients o^, k = 0, • • • ,ft— 1 which will only be satisfied if they vanish. 
In other words, the solution will behave like (x — X+) 71 near the horizon corresponding to 
a pure outgoing wave. However, since ip now vanishes at the horizon, (|2.22j ) implies that 
ijj vanishes everywhere. So there are no nontrivial solutions with iui/K equal to an integer. 
As we saw in section two, this means that if one wanted to include outgoing modes near 
the (future) horizon, the solution would not be smooth there. 

4. Discussion of results 

The numerical procedure described above can be applied to both large black holes 
(r + ^> R) and intermediate size black holes (r+ ~ R). In this section we describe the 
results. We set R = 1, and decompose the quasinormal frequencies into real and imaginary 
parts: 

uj = ur — iuj (4-1) 

With the sign chosen in ( ^4.1| ), uii is positive for all quasinormal frequencies. 

In Table 1., we list the values of the lowest quasinormal mode frequencies for I = and 
selected values of r + , for the four, five, and seven dimensional Schwarzschild-AdS black 
holes. For large black holes, both the real and the imaginary parts of the frequency are 
linear functions of r + . Since the temperature of a large black hole is T = (d— l)r + /4n, it 
follows that they are also linear functions of T. This is clearly shown in fig. 1 and fig. 2, 
where ui and ur respectively are plotted as a function of the temperature for the four, 



11 





4d BH modes 


5d BH modes 


7d BH modes 
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100 


266 3856 


184.9534 


274 6655 


311.9627 


261.2 


500.8 


50 


133.1933 


92.4937 


137.3296 


156.0077 


130.7 


250.4 


10 


26.6418 


18.6070 


27.4457 


31.3699 


26.07 


50.35 


5 


13.3255 


9.4711 


13.6914 


15.9454 


12.96 


25.57 


1 


2.6712 


2.7982 


2.5547 


4.5788 


2.16 


7.27 


0.8 


2.1304 


2.5878 


1.9676 


4.1951 






0.6 


1.5797 


2.4316 


1.3656 


3.8914 






0.4 


1.0064 


2.3629 


0.7462 


3.7174 







Table 1: The lowest quasinormal mode frequency for the 4, 5, and 7 dimensional 
Schwarzschild-AdS black hole for some selected black hole sizes. 




T 



Fig. 1: For large black holes, uij is proportional to the temperature. The 
top line is d = 4, the middle line is d = 5 and the bottom line is d = 7. 



five, and seven dimensional cases. The dots, representing the quasinormal modes, lie on 
straight lines through the origin. In fig. 1, the top line corresponds to the d = 4 case, the 
middle line is the d = 5 case, and the bottom line is the d = 7 case. Explicitly, the lines 
are given by 

wj = 11.16 T for d = 4 
w/ = 8.63T ford = 5 

w/ = 5.47T ford = 7 (4.2) 
Notice from Table 1 that as a function of r+, oji is almost independent of dimension. The 
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Fig. 2: For large black holes, ujr is also proportional to the temperature. 
The top line is now d = 7, the middle line is d = 5 and the bottom line is 
d = 4. 

difference in these slopes is almost entirely due to the dimension dependence of the relation 
between r + and T ( |2.1U[ ). In contrast, ujr does depend on the dimension, and in fig. 2, the 
order of the lines is reversed: 

u R = 10.5 T for d = 7 
ujr = 9.8 T for d = 5 

w fl = 7.75 T for d = 4 (4.3) 

This linear scaling with the temperature is in agreement with the general argument in 
section 2. According to the AdS/CFT correspondence, r = 1/ujj is the timescale for the 
approach to thermal equilibrium. Eq. ( (4.2| ) is one of the main results of this work. 

For the intermediate size black holes, the quasinormal frequencies do not scale with 
the temperature. This is clearly shown in fig. 3 which plots uji as a function of r + for 
d = 4 black holes with r+ ~ 1. To a remarkable acuracy, the points continue to lie along 
a straight line uji = 2.66 r+. The dashed curve represents the continuation of the curve 
ujj = 11.16 T shown in fig. 1. to smaller values of r + . (For large r_|_ these two curves are 
identical.) It is not yet clear what the significance of this linear relation is for the dual 
CFT. Some speculations are given in section 6. Since the quasinormal frequencies can be 
computed to an accuracy much better than the size of the dots in fig. 3, one can check 
that the points actually lie slightly off the line. This is shown more clearly in the five 
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0.5 1 1.5 2 2.5 3 

Fig. 4: u!i for intermediate black holes in five dimensions. The solid line is 
(jJi = 2.75 r + , and the dashed line is ouj = 8.63 T. 

dimensional results in fig. 4. Once again the dashed curve is the continuation of the curve 
ujj = 8.63 T shown in fig. 1, and the solid curve is the line uj = 2.75 r + that it approaches 
asymptotically. 

The real part of the quasinormal frequencies are shown in similar plots in fig. 5 for 
d = 4 and fig. 6 for d = 5. ujr approximates the temperature more closely than the black 
hole size, but it is clear from fig. 6 that it is not diverging for small black holes. 

We have so far discussed only the lowest quasinormal mode with / = 0. We have also 
computed higher modes and modes with nonzero angular momentum, but the numerical 
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0.5 1 1.5 2 2.5 3 

Fig. 6: ujr for intermediate black holes in five dimensions. The solid line is 
ton = 3.12 r + , and the dashed line is ujr = 9.8 T 

accuracy decreases as one increases the mode number n or /. So we restrict our attention 
to relatively small values of n and /. 

For large black holes, in both four and five dimensions, we find that the low lying 
quasinormal modes are approximately evenly spaced in n. In particular, for r+ = 100, 
u>i(n) ~ 41 + 225 n and ivu(n) ~ 54 + 131 n in four dimensions, whereas ivi(n) ~ 73 + 201 n 
and L0n(n) « 106 + 202 n in five dimensions. 

Increasing the angular momentum / mode has the surprising effect of increasing the 
damping time scale ( loi decreases), and decreasing the oscillation time scale {ujr increases). 
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1 2 3 4 5 6 

Fig. 7: Dependence of uj on / for four dimensional black hole with r + 
The smaller points are ur, the larger points are ujj. 



This is shown in fig. 7, where ujr (smaller points) and ujj (larger points) are plotted against 
I for low values of Z.i An important open question is the behavior of uj as I — > oo. It 
appears to decrease with /, but the general argument in section 2 shows that it cannot 
become negative. If ojj continues to decrease with I, then the late time behavior of a 
general perturbation will be dominated by the largest / mode. The large / behavior of 
uji is currently under investigation. Preliminary results indicate that the frequencies stay 
bounded away from zero. If this were not the case, and u>i approached zero fast enough, 
then a general superposition of all spherical harmonics could decay at late times only as a 
power law. However, even this would not be a problem for the AdS/CFT correspondence, 
since the decomposition into spherical harmonics can be done in the boundary field theory 
as well. The statement is that, e.g., a perturbation of < F 2 > with given angular depen- 
dence Yi on S 3 will decay exponentially with a time scale given by the imaginary part of 
the lowest quasinormal mode with that value of I. 



5. Comments on small black holes 

In this section we briefly discuss the extrapolation of the quasinormal frequencies to 
the small black hole regime (r+ <C R). Our numerical approach becomes unreliable in this 
regime, so we cannot compute them directly. Instead, we must rely on indirect arguments. 
But first we give some motivation for exploring this question. 

5 The size of the dots is not related to the accuracy of the calculation. 
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Small AdS black holes are not of direct interest for the AdS/CFT correspondence. This 
is because an extended black hole of the form Schwarzschild-AdS cross S m is unstable 
to forming a black hole localized in all directions whenever the radius of the black hole 
is smaller than the radius of the sphere. This is a classical instability first discussed by 
Gregory and Lafiamme [|16[ . It is quite different from the Hawking-Page transition [|l7l|T8| 
which applies to black holes in contact with a heat bath. In that case, when the black hole 
is of order the AdS radius it undergoes a transition to a thermal gas in AdS. The Hawking- 
Page transition can be avoided if we consider states of fixed energy, not fixed temperature. 
Then black holes dominate the entropy even when r + < R and continue to do so until 



r+/R is less than a negative power of N [l9|,20|j. The situation is very similar to the old 
studies of a black hole in a box. For fixed total energy, the maximum entropy state consists 
of most of the energy in the black hole, and a small amount in radiation. Unfortunately, 
the stable small black hole configuration must be a ten or eleven dimensional black hole 
(by the Gregory-Laflamme instability) and is not known explicitly. 

Nevertheless, there may be other applications of the quasinormal modes of small black 
holes in AdS. One possibility comes from the striking fact (shown in fig. 3) that for d = 4, 
Uj is proportional to r_|_ to high accuracy. As we will discuss below, the slope of this line, 
2.66, turns out to be numerically very close to a special frequency which arises in black 
hole critical phenomena first studied by Choptuik . To explore this possible connection, 
one needs to consider quasinormal modes of small black holes. 

From the intermediate black hole results shown in the previous section, it is tempting 
to speculate that as r + — > 0, u>i — > 0, and u)r — ► constant. Since the decay of the field is due 
to absorption by the black hole, it is intuitively plausible that as the black hole becomes 
arbitrarily small, the field will no longer decay. It is even possible that the quasinormal 
modes approach the usual AdS modes in the limit r + — > 0, although this is not guaranteed 
since the boundary conditions at r = r_|_ do not reduce to regularity at the origin as r_|_ — > 0. 



If they do approach the usual modes in this limit, then our must approach d — 1 pi]] . Of 
course, in the context of string theory, one cannot trust the Schwarzschild-AdS solution 
when the curvature at the horizon becomes larger than the string scale. By taking the AdS 
radius R sufficiently large, one can certainly use this solution to describe some small black 
holes, but the geometry would have to be modified before the limit r + — > is reached. 

It has been shown that the low energy absorption cross section for massless scalars 
incident on a general asymptotically flat spherically symmetric black hole is always equal 
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Fig. 8: The curved line is a fit to the modes of a small black hole in d 
The modes approach the straight line shown at large r+. 



5. 



to the area of the event horizon [[22 . We can use this to estimate the imaginary part of 
the lowest quasinormal mode for a small AdS black hole as follows. Imagine a wave with 
energy of order 1/R propagating toward a black hole with r + <C R. Then the spacetime 
around the black hole is approximately Schwarzschild, and the low energy condition is 
satisfied, so the amplitude of the reflected wave $ r will be reduced from the amplitude 



of the incident wave $i, by 1 — ($ r /$i)" 



rv f 



d-2 
+ ■ 



After a time of order the AdS radius, 



,d-2 
+ ■ 



d-2 



the reflected wave will bounce off the potential at infinity with no change in amplitude. It 
will again encounter the black hole potential and be partly absorbed and partly reflected. 
Repeating this process leads to a gradual decay of the field $ ~ e~ at with a ^ r 

This suggests that for small black holes, u>i should scale like the horizon area r + . 
In the large black hole regime, on the other hand, we know that the modes should scale 
linearly with r + . To check this, we consider a simple ansantz which interpolates between 

ar m 

these two regimes and see how well it fits the data. Consider the function ui^r^) = &+ 
(where a corresponds to the asymptotic slope). For each m we choose b to give the best 
fit to the intermediate black hole data, and see which m yields the lowest overall error (as 
measured by x 2 ). 

In the five dimensional case, using seven points between r + = .4 and r+ = 1, we indeed 
find that m = 3 gives the best fit: x 2 ~ 9 x 10 -6 for m = 3, as opposed to x 2 ~ 3 x 10 -2 
for m = 2 and m = 4. The actual fit, shown in fig. 8 along with the modes and the 
asymptotic line, is given by ujj(r + ) w o 0748+r^" • "'' n ^ our dimensions the story is much less 
clear, since there is no significant difference between the fit with m = 2 and the fit with 
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m = 3. This could be due to the fact that the data for intermediate black holes have not 
yet started to significantly deviate from a straight line. 

To see the possible connection with black hole critical phenomena, consider the evo- 
lution of a self gravitating spherically symmetric scalar field (in an asymptotically flat 4d 
spacetime). It is clear that weak waves will scatter and go off to infinity, just like in flat 
spacetime. Strong waves will collapse and form a black hole. Choptuik [[12 studied one 
parameter families of initial data (labelled by p) which interpolated between these two ex- 
tremes. In each case there is a critical solution p = p* which marks the boundary between 
forming a black hole or not forming one. The late time behavior of this critical solution 
turns out to be universal. It has precisely one unstable mode which grows like e xt with 
A = 2.67. This mode is responsible for the famous scaling of the black hole mass for p just 
above the critical value, Mbh ~ (p — P*) 7 where 7 = 1/A = .374. (For a review, see [23*f.) 

The numerical value of A is very close to the slope, 2.66, of the line in fig. 3 giving the 
imaginary part of the quasinormal mode frequencies. Since both numbers involve imag- 
inary frequencies for spherically symmetric scalar fields in four dimensions, it is natural 
to wonder if there might be a deeper connection between these two phenomena. Unfor- 
tunately, it appears at the moment that the agreement is just a numerical coincidence. 
The first thing one might check is whether the agreement continues in higher dimensions. 
Although the critical solutions for black hole formation have not been studied in five or 
seven dimensions, they have recently been calculated in six dimensions [23] with the result 
A = 1/.424 = 2.36. We have redone our calculation in six dimensions and do not find 
agreement. The slope of uj as a function of r + turns out to be 2.693. Another difference 
is that the exponents in black hole critical phenomena are known to be independent of 
the mass of the scalar field. We have checked that the quasinormal frequencies of large 
and intermediate black holes do depend on the mass. One might expect that if there is a 
connection between these two phenomena, it would apply in the limit of small AdS black 
holes. However, we have seen that the modes of small black holes actually deviate from 
the linear relation, so the significance of the asymptotic slope is not clear. While it is 
still possible that some deeper connection exists (perhaps just in four dimensions and for 
massless fields) it appears unlikely. 

As an aside, we note that if one repeats the calculation of critical phenomena for 
spacetimes which are asymptotically AdS, the late time results will be quite different. 
Since energy cannot be lost to infinity, if one forms a black hole at all, it will eventually 



19 



grow to absorb all the energy of the initial state. 



6. Conclusions 

We have computed the scalar quasinormal modes of Schwarzschild-AdS black holes in 
four, five, and seven dimensions. These modes govern the late time decay of a minimally 
coupled massless scalar field, such as the dilaton. For large black holes, it is easy to see 
that these modes must scale with the black hole temperature T. By the AdS / CFT corre- 
spondence, this decay translates into a timescale for the approach to thermal equilbrium in 
the CFT, for large temperatures and perturbations dual to the scalar field. The timescale 
is simply given by the imaginary part of the lowest quasinormal frequency, r = 
From ( f4.2|) , for perturbations with homogeneous expectation values (/ = modes) these 
timescales are r = .0896/T for the three dimensional CFT, r = .116/T for the four di- 
mensional super Yang-Mills theory, and r = .183/T for the six dimensional (0,2) theory. 
As we mentioned earlier, these time scales are universal in the sense that all scalar fields 
with the same angular dependence will decay at this rate. Perturbations associated with 
other linearized supergravity fields will decay at different rates, given by their quasinormal 
mode frequencies. 

Perhaps the most surprising aspect of our analysis are the results for intermediate size 
black holes. For black holes with size of order the AdS radius, we find that the quasinormal 
frequencies do not continue to scale with temperature, but rather scale approximately 
linearly with horizon radius. We do not fully understand the implications of this linear 
relationship for the dual field theories, but we can make the following comments. If one 
considers the field theory at constant temperature, and slowly lowers the temperature, 
then one encounters the Hawking-Page transition |T^,[r5[. At this point the supergravity 



description changes from the euclidean black hole to a thermal gas in AdS. For these low 
temperatures, the relaxation time might still scale with the temperature, but it cannot 
be computed by a classical supergravity calculation, and is not related to quasinormal 
frequencies. 

To interpret the quasinormal frequencies of intermediate size black holes, we must 
consider a microcanonical description. Consider all states in the CFT with energy equal 
to the supergravity energy. Most of these states will be macroscopically indistinguishable, 
in the sense that they will all have the same expectation values of the operators dual to 
the supergravity fields. If the only nonzero expectation value is the stress energy tensor, 
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the states are described on the supergravity side by just the black hole. If you perturb one 
of these CFT states to one which is macroscopically slightly different, it will decay to a 
typical state with a timescale set by the lowest quasinormal mode. The results in section 
4 show that this decay time is determined by the size of the black hole in the supergravity 
description. Of course the field theory knows about the black hole size since its entropy is 
given by the black hole area. However the fact that, in this range of energy, the frequency 
scales linearly with the radius is puzzling. 

The fact that the quasinormal frequencies do not continue to scale with temperature 
is also interesting for the following reason. For a certain range of energies, the supergravity 
entropy S(E) is dominated by small (ten or eleven dimensional) black holes! [|19|. This 
means that the effective temperature, defined by dS/dE = 1/T, has the property that 
it decreases as the energy increases, i.e. the specific heat is negative. By the AdS/CFT 
correspondence, the same must be true in the dual CFT. (This is not a problem since it 
applies to only a finite range of energies.) If the quasinormal modes continued to scale with 
the temperature, then this negative specific heat would have dynamical effects. Instead, 
we find that the relaxation time increases monotonically with decreasing energy. 
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Appendix: Evaluation of quasinormal modes 

As discussed in section 3, in order to find the quasinormal modes, we need to find 
the zeros of Y^=o a n{.^) (~ x +) n m the complex u plane. We compute the quasinormal 
modes using Mathematica, in the following way. We search for the zeros, u>n, of ip^(uj) = 
^^ =0 a n (c<;) (— x + ) n by looking for the minima of I^at] 2 , and checking that the value at 
the minimum is zero, \i^n(^n) | 2 = 0. (In practice, there are numerical errors in the 
computation, so the value at the minimum is instead ~ 10 -14 , or smaller.) In order to 

6 As we discussed in the previous section, these black holes can be quantum mechanically 
stable, since they are in equilibrium with their Hawking radiation. 
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Fig. 9: Convergence plot for a five dimensional black hole with r + = 0.6 

find the correct minimum, we need to specify an initial guess for ujn- (This sometimes 
poses difficulties in searching for new modes, and apart from using analytical or intuitive 
understanding as our guide, we are forced to resort to trial and error.) How close to the 
actual minimum one is required to start depends on the parameters; for the n = 1, I = 
mode of reasonably-sized black hole, this seldom poses any limitations. 

To obtain an accurate estimate of the quasinormal frequencies u, we typically need to 
compute on the order of N = 100 partial sums, depending on the dimension d, the black 
hole size r+, and the mode (i.e. n and /). Roughly speaking, at a fixed partial sum N, the 
relative error in the computed quasinormal frequency grows as r + decreases, and as /, n, 
or d increases. 

The task of determining the mode to the necessary accuracy is fortuitously simplified 
by the fact that the "convergence curve" has a surprisingly simple form. In particular, 
once the partial sums have converged to sufficient accuracy, the variation of u>n is given by 
an exponentially damped sine as a function of N, i.e. ujn ~ u + ce~ N ^ a sin(WV + d) where 
a and b depend on the physical parameters such as r+, while c and d just depend on which 
partial sum we start with. In fact, we can use a fitting algorithm in Mathematica to fit 
these convergence curves. An example is given in fig. 9, where the dots represent Re{uj^) 
for a particular set of parameters (d = 5, r + = 0.6, n = 1, and / = 0), and the solid curve 
is the corresponding fit. This simplification allows us to determine the mode with a much 
higher accuracy then we would be led to expect from the spread of ujn- It also allows us 
to confirm that the numerical errors in the computation of each ujn are negligible, since 
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otherwise, one would expect a more noisy distribution. 

Often the quickest way to obtain the quasinormal mode is to simply look for the 
minimum of \iPn\ 2 near various initial guesses for the frequency, but when that method 
fails, we can also adopt a more systematic approach; eliminating the possibility of occurence 
of quasinormal modes in a given frequency range. This may be carried out in a more 
systematic manner due to the fact that Re(ipN(w)) and Im(ipN(w)) are conjugate harmonic 
functions of uj, which must satisfy the maximum principle. Thus, if we find that i/jn(oj) is 
bounded inside a given region of the complex uj plane, and either Re(ipN(oj)) or Ii7i(i()n(uj)) 
remains nonzero everywhere on the boundary, then ipw is necessarily nonzero everywhere 
inside that region. This ensures that there can be no quasinormal modes with these 
frequencies. We can thus systematically search for the lowest modes by eliminating the 
low frequency regions until we find the modes. 

Once we find one mode for a given set of parameters, continuity of the solution allows 
us to trace the mode through the parameter space; that is, we can find uj for nearby values 
of r + and /. Also, once we know the n = 1 and n = 2 modes for a fixed r + and /, the equal 
spacing between the modes allows us to find the higher n modes (provided the numerical 
errors stay small). 

Thus, the procedure for finding uj(r + ) is the following: We first consider a large black 
hole, where the convergence is good at a low partial sum, e.g. N = 40. For such a low 
cut-off N on the partial sum, we may easily compute ipN(u},r+,l) in full generality. We 
find the desired mode ujn(t+,ti,1) using the method described above, and we can check 
the convergence by comparing this result with that obtained for the lower partial sums. 
We can now follow the mode to smaller values of r+, until the convergence becomes too 
slow, and we need to compute higher partial sums. It becomes more practical at this point 
to fix all the parameters, and consider iftw as a function of uj only. This has the numerical 
advantage of enabling us to compute the partial sums to much higher N; the drawback, 
of course, is that now we need to recompute the whole series for each r+. 

Note added: 

After this work was submitted, we have extended our computations of the higher / modes for 
large black holes, up to / = 25. From fits of the large I behavior of oji, we find strong evidence 
that the frequencies indeed stay bounded away from zero: In particular, a fit of the form u>i(l) = 
1.12 + - t ^ g has x 2 ~ 2 x 10" 6 , as opposed to \ 2 ~ 10" 3 for a fit with uj t (1 -> 00) -> 0. 
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